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U\ ; ABSTRACT 

^ c*j We calculate multicolor light curves (LCs) of supernovae (SNe) from red supergiants 

p [ (RSGs) exploded within dense circumstellar medium (CSM). Multicolor LCs are cal- 

culated by using a multi-group radiation hydrodynamics code STELLA. If CSM is dense 
enough, the shock breakout signal is delayed and smeared by CSM and kinetic energy 
of SN ejecta is efficiently converted to thermal energy which is eventually released 
\ as radiation. We find that explosions of RSGs arc affected by CSM in early epochs 

when mass-loss rate just before the explosions is higher than ~ 10~ 4 Af Q yr _1 . Their 
characteristic features are that the LC has a luminous round peak followed by a flat 
LC, that multicolor LCs are simultaneously bright in ultraviolet and optical at the 
peak, and that photospheric velocity is very low at these epochs. We calculate LCs for 
various CSM conditions and explosion properties, i.e., mass-loss rates, radii of CSM, 
density slopes of CSM, explosion energies of SN ejecta, and SN progenitors inside, to 
I/*-) ■ see their influence on LCs. We compare our model LCs to those of ultraviolet-bright 

Type IIP SN 2009kf and show that the mass- loss rate of the progenitor of SN 2009kf 
just before the explosion is likely to be higher than 10 -4 M Q yr _1 . Combined with 
the fact that SN 2009kf is likely to be an energetic explosion and has large 56 Ni pro- 
duction, which implies that the progenitor of SN 2009kf is a massive RSG, our results 
indicate that there could be some mechanism to induce extensive mass loss in massive 
RSGs just before their explosions. 
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1 INTRODUCTION 

What kind of stars becomes which kind of supernovae 
(SNe)? There still no conclusive answer to this simple ques- 
tion. Looking into the case of core-collapse SNe (CCSNe) 
which result from death of massive stars, there are several 
discrepancies between answers from observers and theorists. 
Observers get information about the progenitors of CCSNe 
by directly getting images of the SN progenitors by looking 
into the archival images of the SN site |Smarttll2009l . and the 
references therein) or statistically analysing the number of 

* takashi.moriya@ipmu.jp 



SNe appeared (e.g.. [ Smith et all l2010l : iBoissier Sz Prantzosl 
l2009l :[A rcavi et alJl201CT ). On the other hand, theorists have 
been calculating evolution of massive stars for decades to 
get progenitor models of SNe. Theoretical pictures of sin- 
gle star evo l ution without rotation are well summarized in 
iHeeer et all |2003t ). 



One of the big discrepancies raised recently between ob- 
servers and theorists is the zero-age main-sequence (ZAMS) 
mass range of Type IIP SN progenitors. Type IIP SNe are 
the most common type of CCSNe which results from explo- 
sions of red supergiants (RSGs) with a large amount of hy- 
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drogerQ. Both observers and theorists suggest that the mini- 
mum Z AMS mass o f Type IIP SN proge n itors at solar meta l- 
licity is ~ 8 M (jSmartt et all 120091 : iHeeer et all l2003ft . 
Meanwhile, the maximum ZAMS mass of Type IIP SN pro- 
genitors from theoretical calcu lations of solar-me tallicity sin- 
gle star evolution is ~ 25 Mq ijHeeer et al ]|200i and the ref- 
erences therein) but the observa tions suggest that t he max- 
imum ZAMS mass is ~ 17 Mq jSmartt et al.ll2009ft . 

This big discrepancy partly comes from uncertainties in 
theoretical models. One of the big uncertainties in modeling 
stellar evolution is in mass loss. The canonical single star 
modeling uses empirical mass-loss rates but there are many 
mechanisms to induce mass loss which are not taken into 
account in such empirical mass-loss rates. For example, the 
effect of stellar rotations on mass-loss rates are studied for a 
long time and it is found that the maximum ZAMS mass of 
the T ype IIP SNe can be reduced by including stellar rota- 
tions (|Georgy et al.ll2009l . and the references therein). Mass 
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l20ld : IHeeer et al l [i"99~ 

didate for the missing mass-loss mechanism. This kind of 
mass loss due to pulsations occurs in a dynamical timescale 
and it is usually not f ollowed by stellar evolution codes. 
lYoon fc Cantiellol (|2010h analyse the dynamical instabilities 
of RSGs and show that RSGs whose ZAMS mass is above 
~ 17 Mq suffer from this instability and their mass-loss 
rate can potentially be as high as ~ 10 -2 Mq yr _1 . Such 
existence of extensive mass loss in massive RSGs may not 
allow massive RSGs to explode as Type IIP SNe and could 
reduce the maximum ZAMS mass of Type IIP SN progen- 
itors. It is worth noting that the minimum ZAMS mass to 
cause the pulsations is the same as the maximum ZAMS 
mass of Type IIP SN progenitors indicated by observations 
(~ 17 Mq). Nuclear flashes which may occur during the evo- 
lution of massive stars can also be a candidate for the driving 
force of such extensive mass loss i Weaver fe Wooslevl [l979; 
iDessart. Livne. "fc Waldmanllioiol ). If such extensive mass 
loss takes place commonly during the evolution of massive 
stars, theoretical predictions for the maximum ZAMS mass 
of Type IIP SN progenitors can be as low as what observa- 
tions imply. The existence of binary s ystems is also sugg ested 
to be account for this problem (e.g.. ISmith et al1l201C ). 

There are also observed potential Type IIP SN progen- 
itors (RSGs) which are losing their mass with very high 
mass-loss rates. For example, a RSG VY Canis Majoris 
is estimated to be losing its mass with the rate 1 — 2 x 
10 -3 Mq yr _1 from the dir ect observations of circumste l- 
lar medium (CSM) around it (|Smith. Hinkle fc Rvdel200gft . 
Observations of dusts around another RSG IRAS05280-6910 
also indicate the extensive mass loss (~ 1(T 3 Mq yr- 1 ) of 
the RSG jBover et al.ll2010ft . Some Type IIP SNe are re- 
ported to show the possible effect of dense CSM in their 
LCs and spectra. A representative exampl e of this kind is 
ultra violet (UV)-bright Type IIP SN 2009kf (|Botticella et all 
l2010ft . SN 2009kf was bright in UV for ~ 10 days during 
its early epochs as well as in optical. Later, the LC trans- 
formed to that of a Type IIP SN and the spectra taken 
at later epochs are classified as Type IIP. The fact that SN 
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2009kf was bright in UV and optical at the same time makes 
it difficult to explain SN 2009kf without CSM interaction. 
This is because usual Type IIP SNe cool adiabatically in 
their early epochs and they become optically bright after 
UV b rightness gets weaker. lUtrobin. Chugai. fc Botticellal 
(|201(]ft show a LC model of SN 2009kf without CSM in- 
teraction but the required explosion energy is found to be 
very high (2.2 x 10 52 erg). Another Type IIP SN 2007od 
l| Andrews et al.l l2010ft had late phase spectra similar to 
those of Type Iln SN 1 998S, which are c onsidered to result 
from CSM interaction (|Pozzo et al.ll2004h . Also, SN 2007od 
showed the possibility of the existence of light echos in its 
LC and it may indicate the exi stence of CSM around the 
progenitor jjAndrews et al.ll2010ft . 

The existence of dense CSM around SN progenitors 
at the time of their explosions can affect their LCs. If a 
SN explosion occur in dense CSM, the shock breakout sig- 
nal from the SN diffuse out in CSM. In addition, the SN 
ejecta is decelerated by the dense CSM. Thus, the kinetic 
energy of the ejecta is efficiently converted to the ther- 
mal energy which is eventually released as radiation en- 
ergy. This extra heat source makes the SN brighter than 
usual SNe. So far, SNe powered by the interactions of SN 
ejecta and dense CSM are mainly modeled in the context 
of Type Iln SNe showing narrow emission lines which are 
a cle ar evidence of the interaction of SN ejecta and CSM 
~ l[200l;IWooslev. Blinnikov. fc Hegeij|2007l ; 
Hvan Marie et all l201fj| ; iMetzgeij I201C 



(e.g..lChugai et al 



IDessart et ail 2009 



Here, given the possibilities that RSGs can have enhanced 
mass-loss rates and that dense CSM can exist around RSGs 
at the time of the explosions, we investigate the effect of 
dense CSM on L Cs from SNe of RSGs wit h dense CSM. 
Some works (e.g., iFalk fc ArnettJll973l , Il977f ) have already 
done with the similar conditions but, in this paper, we 
model them more systematically and with better treatments 
of physics. To model the LCs powered by the shock in- 
teractions, we have to follow radiation hydrodynamics be- 
cause the energy source of the radiation is a hydrodynamical 
shock wave. Thus, we adopt a multi- group radiation hydro- 
dynam i cs code STELLA developed by iBlinnikov fc Bartunovl 
l| 1993ft : IBlinnikov etafl (|l99a 12006ft ■ which has been ap- 
plied for modeling LCs powered by shock interactions 
fe.gJChugai et alj|2004l;IWooslev. Blinnikov. fc Hegerll2007l : 
IBlinnikov fc Sorokinall2010ft . 

This paper is composed of five sections. We first briefly 
summarize STELLA and the pre-SN models we adopt in Sec- 
tion [2] Synthesized LCs with various CSM conditions are 
shown in Section In Section 21 we compare our models 
with the LCs of SN 2009kf and show that the LCs powered 
by the interaction of SN ejecta and dense CSM is actually 
consistent with the LCs of SN 2009kf. We also get some im- 
plications for the CSM around the progenitor of SN 2009kf 
and what kind of mass loss should have happened before the 
explosion of SN 2009kf. Conclusions are given in Section [S] 



2 STELLA CODE AND PRE-SUPERNOVA 
MODELS 

STELLA is a on e- dimensional multi-group radiation hydrody- 
nami c s cod e (|Blinnikov fc Bartunovl [19931 : IBlinnikov et al.l 
1 19981 . 2006). STELLA calculates the spectral energy distri- 



© 2011 RAS, MNRAS 000,mH6l 



Supernovae from RSGs with Extensive Mass Loss 3 



Table 1. List of LC models without CSM 



Name 


Progenitor 


Explosion Energy 
10 51 erg 


Radiation Energy 
10 49 erg 


sl3e3 


sl3 


3 


3.0 


sl5el 


slS 


1 


1.3 


sl5e3 


sl5 


3 


3.2 


sl5e5 


sl5 


5 


4.9 


sl8e3 


sl8 


3 


3.7 


s20e3 


s20 


3 


1.2 



Radiation energy emitted in 50 days since the explosion. 



butions (SEDs) at each time step and we can get multi- 
color LCs by convolving filter functions to the SEDs. All 
the calculations are performed by adopting 100 frequency 
bins from 1 A to 5 x 10 4 A. STELLA implicitly treats time- 
dependent equations of the angular moments of intensity 
averaged over a frequency bin. Local thermodynamic equi- 
librium is assumed to determine the ionization levels of ma - 
terials. For the details of STELLA, see lBlinnikov et al.1 (|200fj| ') 
and the references therein. Comparisons of STELLA with 
other numer ic al codes are pr o vided in, e.g.. iBlinnikov et al.l 
GUI, 120031 ); IWooslev et all l|2007fl and analytical models 
are also compared to the n umerical results of STELLA (e.g., 
iRabinak fc Waxmanl [201ll ). In this paper, explosions are 
treated as a thermal bomb by injecting thermal energy just 
below the mass cut 1.4 Mq. We do not follow explosive nu- 
cleosynthesis because we focus on the early epochs when the 
effect of explosive nucleosynthesis on SN LCs is negligible. 
We construct pre-SN models by attaching 

CSM to the progen itor models calculated by 

IWooslev. Heger, fc Weaver] ||2002|). As th e progenitor 
models of IWooslev. Heger. fc Weaver! (2002) do not take 
into account of the extensive mass loss we are interested 
in, we artificially attach CSM to the outer most layer of 
the progenitor models. The composition of CSM is solar 
metallicity and is the same as the surface of the RSG 
models which we adopt. Compared to the effect of the 
CSM parameters, the LCs are less affected by the RSG 
models adopted inside (Section I3.2.6[) . This justifies our 
simple way to construc t the pre-SN models. Among th e 
pre-SN models shown in lWooslev. Heger. fc Weaver! (|2002h . 
we use the solar- metallicity single star models of RSGs; 
sl3, sl5, sl8, and s20. The ZAMS masses of the models 
sl3, sl5, sl8, and s20 are 13 M©, 15 M , 18 M Q , and 20 
Mq, respectively. If they are exploded without CSM, their 
LCs show a long plateau phase and thus they are Type IIP 
SN progenitors (Section 13. ip . Although we do not follow 
explosive nucleosynthesis, 56 Ni is in the core of the pre-SN 
models due to the nuclear statistical equilibrium established 
in the core. However, at the early epochs we are interested 
in (first ~ 50 days since the explosions), the photons 
originated from 56 Ni decay do not leak from the ejecta so 
m uch and contribute little to LCs. According to Figure 2 
of lKasen fc Wooslevl (|2009l ). the effect of 56 Ni on Type IIP 
SN LCs appear after ~ 50 days since the explosions and 
thus our LC models are not applicable from around that 
epoch. Density structures of the pre-SN models are shown 
accordingly in the following sections. 




20 40 60 80 100 

days since the explosion (rest frame) 

Figure 1. Bolomctric LCs of the SNe from the progenitors with- 
out CSM. 



3 LIGHT CURVES 

In this section, we present the LCs calculated by STELLA. 
First, we show LCs of the explosions of the progenitors with- 
out CSM for references (Section I3.1|l and then LCs of the 
explosions with CSM (Section [ 



3.1 Explosions without CSM 

In this section, we present the LCs for the progenitors with- 
out CSM. The aim of these calculations is to provide ref- 
erences to be compared with LCs with CSM. Note that, as 
mentioned in Section [2] the LCs shown here are not appli- 
cable after ~ 50 days since the explosions because we do not 
follow explosive nucleosynthesis. Table[T]is the list of the LCs 
shown in this section. Previous studies also calculated the 
multicolor LCs of Type IIP SNe from th e same progenitor 
mode l s with different numerical codes jKasen fc Wooslev 
| 2009l ; iDessart. Livne. fc Waldmanl l20ld ; IDessart fc Hillier 
120111 ) and from di fferent progenitor models with th e same 
numerical code (iBaklanov. Blinnikov. fc Pavlvukl 120051 ; 
iTominaga et al.ll2009l . l201lh . 

In Figure [T] the bolometric LCs with the same explosion 
energy (3 x 10 51 erg) are presented. When a shock emerges 
from a stellar surface, a bolometric LC is suddenly bright- 
ened due to the shock breakout. Then, the ejecta cools adi- 
abatically. When the outer layer of the ejecta reaches the 
recombination temperature of hydrogen, the LC becomes 
flat until the photosphere reaches the bottom of the hydro- 
gen layer (plateau phase). After the plateau phase, the LC 
follows the decay line of 56 Co which existed as 56 Ni in the 
core of the pre-SN model even though we do not calculate 
the explosive nucleosynthesis (Section Note that our LC 
models are not applicable after the late epochs of the plateau 
phase. 

The bolometric LCs with the same progenitor (sl5) but 
different explosion energies are shown in the top left panel 
of Figure [2] In the other panel of Figure [2] we also show the 
multicolor LCs of each model. Optic al LCs are ob tained by 
applying the Bessell UBVRI filters <|Bessel]||l99oT ) and UV 
LCs are derived by using the near-UV (NUV) imaging filter 
of the Galaxy Evolution Explorer (GALEX) satellite whose 
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days since the explosion (rest frame) 




40 60 80 100 120 
days since the explosion (rest frame) 



160 




40 60 80 100 120 140 
days since the explosion (rest frame) 



40 60 80 100 120 
days since the explosion (rest frame) 



160 




160 



Figure 2. Bolometric (top left) and multicolor (the others) LCs of the SNe from the progenitor sl5 without CSM. Models with several 
explosion energies are shown. 



central wavelength is around 2300 A l]Morrissev et al.ll2005l ). 
After the first brightening due to shock breakout, NUV first 
become bright because of the adiabatic cooling of the ejecta. 
Then, LCs become bright in optical as NUV become fainter. 
This is an important feature of usual Type IIP SNe. They are 
not bright in UV and optical at the same time. The behavior 
of the multicolor LCs of the other progenitors is qualitatively 
the same as the multicolor LCs shown in Figure [2] 

3.2 Explosions with CSM 

In this section, we investigate the effect of CSM around 
RSGs on SN LCs. After the explosion of the progenitor, 
a shock wave propagates inside the progenitor. At these 
epochs, CSM is not ionized and optically thin. Then, after 
the shock wave has gone through the surface of the pro- 
genitor, a precursor wave appears in the CSM and CSM is 
ionized. The precursor wave propagates very fast and soon 
reaches to the surface of CSM. The optical depth (tcsm) of 
CSM after the shock wave reached the surface of the pro- 
genitor is expressed as 

t"csm = / K(r)p(r)dr~ / n(r)p{Ro) (jr) dr > (!) 

J R J R ...in, 

where k is opacity, p is density, R is the CSM radius, Ro 
is the radius of the inner edge of CSM where CSM is con- 



nected to the progenitor inside, i.e., the radius of the pro- 
genitor inside, and Ri is the ionization front in CSM. tcsm 
depends on the thickness of the ionized layer because « 
above the ionization front is very low (k ~ 10 -4 cm 2 g _1 ). 
In order to estimate a condition in which CSM can af- 
fect LCs, we assume that CSM below Ri is fully ionized 
and the Thomson scattering is the predominant source of 
opacity (ft(r) = 0.33 cm 2 g _1 with solar metallicity). With 
Ro ~ 5 X 10 13 cm (Figure [3] Section GLMJ , Ro < Ri, and 
a > 1, the condition to be tcsm > 1 is 

p(R ) > 6x HT 14 (a-l) g cm" 3 . (2) 

Thus, when the density of CSM at the radius where CSM 
is connected to the progenitor model inside is more than 
~ 10 -13 g cm -3 , CSM becomes optically thick and LCs are 
expected to be affected by CSM. CSM with mass-loss rates 
higher than ~ 10~ 4 Mq yr _1 satisfies this condition, assum- 
ing that the CSM velocity is 10 6 cm s" 1 (Figure 

In the following sections, we also investigate the depen- 
dence of LCs on several physical parameters of CSM and 
progenitors. The parameters of CSM adopted are mass-loss 
rates [10~ 4 - 10~ 4 M© yr" 1 ], radii of the outer edge of CSM 
(CSM radii) [5 x 10 14 - 3 x 10 15 cm], and density slopes 
(p oc r~ a ) [a = 2,1.5]. Several progenitor models inside 
are also adopted [sl3, sl5, sl8, s20]. The density structures 
of the pre-SN models with CSM are shown accordingly in 
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Figure 3. Density structures of the pre-SN models with different 
mass-loss rates and density slopes. The density structures shown 
with solid lines have CSM density slope of p oc r~ 2 and those 
shown with dashed lines have CSM density slope of p oc r~ 1 . 
The pre-SN models are constructed by attaching CSM to the 
progenitor model sl5. 




10 20 30 40 50 60 
days since the explosion (rest frame) 

Figure 4. Bolometric LCs of the models with different mass- 
loss rates. The models shown are sl5w2r20mle3 (10 -1 Mq yr - 1 ), 
Sl5w2r20m2e3 (10~ 2 Mq yr" 1 ), sl5w2r20m3e3 (1CT 3 Mq yr" 1 ), 
and Sf5w2r20m4e3 (10 -4 Mq yr -1 ). The bolometric LC of the 
model s!5e3 (no CSM) is also shown for comparison. 



the following sections. The model parameters and results are 
summarized in Tables [2] and [3] 



3.2.1 Effect of CSM 

Typical bolometric LCs which are affected by CSM are 
shown in Figure [4] with the bolometric LC from the model 
without CSM (sl5e3). We focus on the model sl5w2r20m3e3 
to describe the effect of CSM in this section. CSM mainly 
affects the LCs at the early epochs, roughly until the sud- 
den drop in the LCs which can be seen at around 25 days 
in the model sl5w2r20m3e3 (1CT 3 Mq yr -1 ). At first, LCs 
have round shapes (until ~ 15 days in sl5w2r20m3e3). The 
round phase is followed by a fiat LC which lasts until the 
sudden drop (~ 15 — 25 days in sl5w2r20m3e3). LCs are 
mainly powered by the interaction of SN ejecta and CSM at 
these epochs. We call these epochs as an interaction-powered 
phase (IPP). 

The IPP appears in the LCs of the SNe with dense 
CSM. When the temperature and CSM density are high 
enough, CSM becomes optically thick and the photo- 
sphere locates in CSM. Photosphere during the IPP is in 
CSM. This can clearly seen in Figure [5] Looking into the 
model sl5w2r20m3e3, the photospheric velocity is at first 
10 6 cm s _1 , which is the initial CSM velocity. This is a char- 
acteristic feature of the explosions within dense CSM. Then, 
the photospheric velocity increases due to the acceleration 
of CSM by the precursor wave. At around ~ 15 days in 
sl5w2r20m3e3, the photosphere reaches the shell between 
SN ejecta and CSM (Figure [SJ. Then, the photosphere goes 
into SN ejecta at ~ 25 days. 

All the LCs affected by CSM have round shapes at first 
(Figure [4j. We briefly discuss why the round phase appears. 
As our models have tqsm > 1, photons cannot escape freely 
from CSM. In addition, our models satisfy the following con- 
dition when the shock wave is propagating in CSM: 




10 20 30 40 50 60 70 80 
days since the explosion (rest frame) 



Figure 5. Evolutions of photospheric velocity of the models 
sl5w2r20m2e3, sl5w2r20m3e3, and sl5e3. The photosphere is de- 
fined as the location where the Rosseland optical depth becomes 
2/3. Observed Ho line velocity of SN 2009kf is also plotted in the 
figure. The explosion date of SN 2009kf is set as the same as in 
Figure [TdI 

f* c 
t s = K(r)p(r)dr < — , (3) 

Jr. v ° 

where R 3 is the radius of the shock wave and v s is the shock 

velocity. Therefore, photons can diffuse out from the shock 

wave and a precursor wave propagates ahead of the shock 

wave (see the left panels of Figure [6j . As the shock velocity 

is typically ~ 10 9 cm s -1 when the shock wave reach Ro, 

c/v s is typically ~ 10 at that time. Thus, photons in the 

models with 

p(R ) < 6 x 1CT 13 g cm" 3 (4) 

start to leak photons from the shock wave just after the 
shock wave reach Tin- The other models satisfy Equation 
when the shock wave propagates in CSM because the de- 
celeration of the shock wave makes c/v s higher and the 
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Table 2. List of LC models with CSM 



Ncii lie 


Progenitor 


Explosion Energy 


Mass-Loss Rate a 


Radius 


Density Slope 


CSM Mass 


Radiation Energy^ 






10 51 ere 


Mr\ vr — 1 


10 15 cm 


ex. oc v a ) 


Mr; 


in 49 erg- 


Sl3w2r20m2e3 


sl3 


3 


10~ 2 


2 


2 


0.65 


20 


Sl5w2r05m2e3 


sl5 


3 


10- 2 


0.5 


2 


0.15 


9.8 


sl5w2r05m3e3 


sl5 


3 


10" 3 


0.5 


2 


0.015 


3.8 


Sl5w2rl0m2e3 


si 5 


3 


1(T 2 


1 


2 


0.31 


13 


sl5w2rl0m3e3 


sl5 


3 


10- 3 


1 


2 


0.031 


4.2 


sl5w2r20m2el 


sl5 


1 


io- 2 


2 


2 


0.65 


5.6 


Sl5w2r20mle3 


sl5 


3 


io- 1 


2 


2 


6.5 


75 


sl5wl.5r20mle3 


sl5 


3 


- 


2 


1.5 


6.5 


80 


Sl5w2r20m2e3 


sl5 


3 


IO" 2 


2 


2 


0.65 


19 


sl5wl.5r20m2e3 


sl5 


3 


- 


2 


1.5 


0.65 


15 


Sl5w2r20m2e5 


sl5 


5 




2 


2 


0.65 


32 


Sl5w2r20m2e7 


sl5 


7 


IO" 2 


2 


2 


0.65 


45 


slowzrzUmoeo 


sl5 


3 


1 n— 3 


2 


2 


0.065 


5.0 


sl5wl.5r20m3e3 


sl5 


3 




2 


1.5 


0.065 


4.8 


Sl5w2r20m4e3 


sl5 


3 


10- 4 


2 


2 


0.0065 


3.4 


Sl5wl.5r20m4e3 


sl5 


3 




2 


1.5 


0.0065 


3.4 


Sl5w2r30m2e3 


sl5 


3 


IO- 2 


3 


2 


0.98 


23 


sl5w2r30m3e3 


sl5 


3 


IO" 3 


3 


2 


0.098 


5.6 


Sl8w2r20m2e3 


sl8 


3 


IO" 2 


2 


2 


0.65 


17 


S20w2r20m2e3 


s20 


3 


IO" 2 


2 


2 


0.65 


17 



a Derived by assuming a constant CSM velocity 10 6 cm s 1 . 
^ Radiation energy emitted in 50 days since the explosion. 



Table 3. Properties of SNe from sl5 at some selected epochs 



iboll (IO 43 erg s- 1 ) v p £ (km s" 1 ) T BE £ (IO 3 K) 

Epoch£ (days) 10 20 30 40 10 20 30 40 10 20 30 40 



sl5el 


0.27 


0.26 


0.23 


0.23 


6900 


6500 


5900 


5000 


15 


9.8 


7.9 


7.0 


sl5e3 


0.84 


0.65 


0.60 


0.57 


12000 


11000 


9500 


7800 


16 


9.8 


8.0 


7.1 


sl5e5 


1.3 


1.0 


0.90 


0.81 


15000 


13000 


12000 


9200 


15 


9.9 


8.0 


7.1 


sl5e7 


1.7 


1.1 


1.2 


1.0 


17000 


16000 


13000 


10000 


16 


9.9 


8.1 


7.1 


Sl5w2r05m2e3 


1.1 


0.63 


0.73 


0.57 


8100 


8300 


8300 


7800 


15 


9.0 


7.7 


7.2 


sl5w2r05m3e3 


0.82 


0.69 


0.61 


0.56 


10000 


10000 


9400 


7800 


13 


9.6 


8.0 


7.2 


Sl5w2rl0m2e3 


22 


0.78 


0.73 


0.64 


950 


7800 


7500 


7600 


42 


9.1 


7.9 


7.2 


Sl5w2rl0m3e3 


1.5 


0.72 


0.61 


0.56 


9800 


9800 


9400 


7800 


11 


9.1 


8.0 


7.2 


Sl5w2r20m2el 


0.81 


3.2 


1.2 


0.67 


14 


86 


130 


150 


15 


21 


14 


10 


Sl5w2r20mle3 


0.049 


1.9 


47 


36 


10 


15 


410 


1200 


6.0 


13 


37 


32 


sl5wl.5r20mle3 


0.064 


1.5 


45 


46 


10 


12 


400 


1400 


6.0 


12 


36 


34 


Sl5w2r20m2e3 


9.1 


8.4 


2.5 


0.92 


65 


500 


7200 


6800 


29 


23 


11 


7.2 


sl5wl.5r20m2e3 


5.0 


7.9 


3.0 


0.64 


43 


290 


6300 


5600 


23 


2.5 


16 


7.5 


Sl5w2r20m2e5 


29 


7.6 


1.2 


1.1 


670 


1000 


8900 


9000 


39 


20 


7.9 


7.3 


Sl5w2r20m2e7 


58 


8.1 


1.7 


1.3 


570 


11000 


11000 


10000 


42 


16 


8.2 


7.4 


Sl5w2r20m3e3 


2.1 


1.3 


0.62 


0.57 


420 


9400 


9200 


7800 


18 


9.5 


8.0 


7.2 


Sl5wl.5r20m3e3 


2.1 


1.5 


0.62 


0.36 


270 


9200 


9200 


7800 


18 


9.8 


8.1 


7.2 


Sl5w2r20m4e3 


0.96 


0.74 


0.61 


0.57 


11000 


10000 


9400 


7800 


16 


10 


8.1 


7.1 


Sl5wl.5r20m4e3 


0.90 


0.77 


0.61 


0.56 


11000 


10000 


9500 


7800 


17 


11 


8.1 


7.2 


Sl5w2r30m2e3 


2.8 


15 


4.4 


2.5 


21 


210 


300 


330 


20 


26 


16 


10 


sl5w2r30m3e3 


3.5 


1.3 


1.1 


0.61 


30 


9400 


9100 


8000 


21 


9.8 


8.7 


7.3 



Bolomctric luminosity. 

Photospheric velocity. Photosphere is where the Rosseland optical depth is 2/3. 

Color temperature at photosphere. Photosphere is where the Rosseland optical depth is 2/3. 

Epoch since the explosion. 
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Figure 6. Physical structures of the models sl5w2r20m2c3 and sl5w2r20m3e3 at two epochs. The epochs are the days since the 
explosion. Black lines show the density structure (left i/-axis). Blue lines are the velocity scaled by 10 s cm s _1 (right j/-axis), purple 
lines are the logarithm of the absolute value of luminosity scaled by 10 40 erg s — 1 (right y-axis), green lines are the logarithm of the 
temperature in Kelvin (right j/-axis), and red lines are Rosseland optical depth measured from the outside (right j/-axis). 



propagation of the shock wave reduces t s . Once Equation 
[3] is satisfied, photons start to leak out from the shock wave 
and this phenomenon is usually observed as shock break- 
out. However, since there is dense CSM at the time of shock 
breakout, photons diffuse out in CSM and shock breakout 
signal become longer compared with the explosions without 
CSM (see also lFalk fc ArnettJll973l . Il977r i. 

In addition, SN ejecta is decelerated by the dense CSM. 
Dense CSM is massive and it has the density structure 
p oc r~ 2 if it is from steady mass loss. Thus, the shock 
wave between the SN ejecta and CSM is decelerated and 



kinetic energy of the SN ejecta is converted to thermal en- 
ergy which is released as radiation energy. As a result, SNe 
with dense CSM emit more photons and become brighter 
tha n those without dense CS M. For further discussion, see, 
e.g.. IChevalier fc Irwinl l|201ll ). Comparing radiation energy 
emitted during early epochs (Tables [T] and [2}, it is clear that 
the effect of the SN ejecta deceleration is dominant radiation 
source during the IPP. In other words, the round phase is 
not just due to the elongation of the shock breakout signal 
seen in the models without CSM (e.g. sl5e3). 

The difference in the rising times and durations of the 
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10 20 30 40 50 60 70 10 20 30 40 50 60 70 

days since the explosion (rest frame) days since the explosion (rest frame) 

Figure 7. Multicolor LCs of the models with different mass-loss rates (solid lines). The dashed lines are the multicolor LCs of the model 
sl5e3 for comparison. LCs of the same color are from the same filter. Note that top left panel has different y-axis range from the other 
panels. 



round LCs during the IPP comes from the difference in dif- 
fusion timescales of CSM (Figure[4]). The models with higher 
mass-loss rates have denser CSM and thus longer diffusion 
timescales. This difference due to diffusion timescales can 
be clearly seen in Figure [6] The upper panels and the lower 
panels represent the same epoch of the explosion with differ- 
ent CSM densities, i.e., different CSM diffusion timescales. 
As is indicated by the temperature waves in CSM which 
are pushed by the photon diffusion, photons in CSM with 
shorter diffusion timescales diffuse out more quickly into 
CSM. Therefore, the rising times and the durations of the 
round phase in the IPP are shorter for the models with the 
smaller CS M diffusion time s cales. Those differences are also 
discussed in lFalk fc ArnettJ (|l973l . I1977T I. 

The round LC in the IPP is followed by the flat LC 
which lasts until the sudden drop of the LC (between around 
15 days and 25 days in sl5w2r20m3e3). During this flat 
phase, the photosphere locates at the dense shell between SN 
ejecta and CSM. The photospheric velocity does not change 
so much during the flat phase (Figure [5]). After CSM above 
the dense shell have become optically thin, the photosphere 
remains at the dense shell until the temperature and density 
of the shell becomes low enough to be optically thin. 

One of the clear characteristics of the LCs with CSM 
is the sudden drop in the LCs. The time of this sudden 



drop corresponds to the time when the dense shell becomes 
optically thin and the photosphere proceeds inward to SN 
ejecta. This can also be seen in the photospheric velocity 
evolution (Figure [5|. The brightness can drop as low as the 
LC without CSM because now the photons are emitted from 
SN ejecta and the physical conditions are the same as those 
of SNe without dense CSM. However, the brightness is still 
slightly more luminous than the LC without CSM for several 
days after the sudden drop. This could be because of the 
extra-heating due to shock and/or deceleration of SN ejecta 
by CSM which makes the adiabatic cooling of SN ejecta less 
efficient. 

There are many differences between our models (RSG 
+ CSM) and RSGs with extended envelopes. First of all, 
it is difficult to have RSGs ext ended to ~ 10 15 cm (e.g., 
IWooslev. Heger. fc Weaverl2002h . What is more, density and 
temperature in RSGs are much higher in RSG envelopes 
than in CSM. This does not allow the shock wave in RSGs 
to satisfy Equation [3] until the shock wave reaches the sur- 
face of RSGs. In other words, the shock wave does not break 
out until it reaches near the surface. On the contrary, in our 
models, Equation[3]is satisfied inside CSM and the precursor 
wave propagates ahead of the shock wave. In addition, the 
kinetic energy of the SN ejecta efficiently converted to ra- 
diation energy because of the deceleration of the SN ejecta 
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Figure 8. Density structures of the pre-SN models with dif- 
ferent radii. Models with solid lines have the mass-loss rate 
10~ 2 Mq yr" 1 and those with dashed lines have 10~ 3 Mq yr" 1 . 
The pre-SN models are constructed by attaching CSM to the pro- 
genitor model s!5. 



by the dense CSM. This precursor wave due to the shock 
breakout in CSM and an additional heating source causes 
the IPP phase. 



3.2.2 Dependence on Mass-Loss Rate 

In this section, we show the effect of mass-loss rates on LCs. 
We adopt the mass-loss rates of 1CP 1 , 10~ 2 , 10 -3 , and 10 -4 
Mq yr~ (the corresponding model names are sl5w2r20ml, 
sl5w2r20m2, sl5w2r20m3, and sl5w2r20m4, respectively). 
For the case of KT 1 M Q yr" 1 , CSM mass is 6.5 Mq and the 
sum of the mass of RSG inside and CSM exceeds the Z AMS 
mass of the progenitor. This is unrealistic but we show the 
results just to see the effect of CSM. Every CSM in the 
models is optically thick (Equation[2] Figure[3]). These mass- 
loss rates are derived by assuming that the CSM velocity is 
10 6 cm s _1 . However, the escape velocity of the surface of 
sl5 is 7.6 x 10 6 cm s -1 . This means that the CSM velocity 
can be higher than 10 cm s~ x at least at the late stages of 
the evolution of the progenitor and thus the actual mass- 
loss rate for a given CSM density could be higher than the 
values we show. The flow from the progenitor may not be 
even steady. However, we assume that CSM results from the 
steady flow from the progenitor with the velocity 10 6 cm s _1 
for simplicity because the CSM velocity has little effect on 
LCs. To see the effect of the mass-loss rates on the LCs, we 
fix the radius of CSM to 2 x 10 cm and the density slope to 
p oc r~ 2 in this section. Also, the explosion energy and the 
progenitor of the SNe are fixed to 3 x 10 51 erg and sl5. The 
density structures of the progenitors with CSM are shown 
in Figure [3] 

Figure [4] shows the bolometric LCs of the SNe with dif- 
ferent mass-loss rates. Since the diffusion timescale of CSM 
becomes larger for denser CSM, rising times and durations of 
round phases in LCs are longer for the models with higher 
mass-loss rates. The maximum luminosity of the LCs be- 
comes larger as CSM becomes denser. This is because the 
shock wave is more decelerated by CSM due to the more 
massive CSM, i.e., more kinetic energy is converted to the 
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Figure 9. Bolometric LCs of the models with different CSM 
radii. The top panel shows the LCs from the models with the 
mass-loss rate 10 -2 Mq yr _1 and the bottom panel shows those 
with 10 — 3 Mq yr — x . The bolometric LC of the model sl5e3 is 
also shown for comparison. 



thermal energy and thus radiation energy. The radiation en- 
ergies of the models are summarized in Table [2] 

The multicolor LCs of the models in this section are 
shown in Figure [7] Each LC is plotted with the multicolor 
LCs of no CSM model sl5e3. It is clear that the LCs dur- 
ing the IPP become bright especially in short wavelengths 
because of high temperatures at the photosphere. 



3.2.3 Dependence on Radius 

In this section, the effect of the CSM radius on LCs is 
investigated. To see the effect of the CSM radius, we fix 
the following parameters; the explosion energy (3 x 10 51 
erg), the density slope (a = 2), and the progenitor in- 
side (sl5). We try two mass-loss rates, 10 -2 Mq yr _1 and 
10 -3 Mq yr -1 . We adopt four CSM radii to see the effect, 
i.e., 5 x 10 14 cm (sl5w2r05m2 and sl5w2r05m3), 1 x 10 15 cm 
(sl5w2rl0m2 and sl5w2rl0m3), 2 x 10 
and sl5w2r20m3), and 3 x 10 15 
Sl5w2r30m3). With the constant CSM velocity 10 6 cm s" 1 , 
the mass loss in the models lasts 16 years, 32 years, 64 years, 
and 96 years, respectively. The density structures of the pre- 
SN models are shown in Figure [8] 

Figure [9] shows the bolometric LCs of the models with 



15 cm (sl5w2r20m2 
cm (sl5w2r30m2 and 
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Figure 10. Multicolor LCs of the models with different radii (solid lines). The dashed lines are the multicolor LCs of the model sl5e3 
for comparison. LCs of the same color are from the same filter. Note that the scale of y-axis in the right bottom panel (sl5w2r30m2e3) 
is different from those in the other panels. 



different CSM radii. The durations of the round phases in 
the IPP are longer for the models with the larger CSM ra- 
dius. This is because the diffusion times of the LCs are longer 
for the models with the larger CSM radius. On the other 
hand, the maximum luminosities of the round phases de- 
crease as the CSM radius increases. As the explosion energy 
and the density structure are similar in each model, the ra- 
diation energy released by the shock interaction is also close 
to each other (see Table [2] for the radiation energy emit- 
ted in each model). Therefore, the difference in the maxi- 
mum luminosities is caused by the difference in the diffu- 
sion timescales. Even if the same energy is released in the 
same timescale, photons are more scattered and distribute 
more uniformly in CSM for models with longer diffusion 
timescales. Thus, the luminosity, i.e., radiation energy re- 
leased from the CSM surface in a unit time, becomes lower 
for the models with larger CSM radii. The flat phase of IPP 
is also longer for the models with the larger radius. 

The multicolor LCs of the models with 10 -2 Mq yr _1 
are shown in Figure 1101 Although sl5w2r05m2e3 has the 
brightest peak bolometric luminosity among the models, the 
NUV and UBVRI band LCs of the model are the faintest. 
This is because the more compact CSM is, the hotter CSM 
becomes. The photosphere of the model sl5w2r05m2e3 is too 
hot during the IPP to emit the radiation in the NUV and 



UBVRI bands. The SEDs of the models at the bolometric 
peak are shown in Figure [TT] 



3.2.4 Dependence on Density Slope 

The CSM density slope of p oc r~ 2 results from the steady 
flow from the central progenitor. However, if mass loss is 
not steady, the density slope does not necessarily have the 
density structure p oc r~ 2 and can be shallower (a < 2) or 
steeper (a > 2). For example, the shallower density slope 
of CSM is sug gested in the mod eling of the spectra of Type 
Iln SNe (e.g JChugai et alj2004 ). We calculate the LCs with 
the slope a = 1.5. As we are fixing the CSM velocity, mass- 
loss rates should change with time to have the density slope 
a = 1.5. The CSM radius (2 x 10 15 cm), the explosion en- 
ergy (3 x 10 51 erg), and the progenitor inside (sl5) are fixed 
in this section. To see the effect of density slopes, we calcu- 
late the models with the same CSM masses as the models 
with the density slope a = 2, i.e., 6.5 Mq (sl5wl.5r20ml), 
0.65 M Q (sl5wl.5r20m2), 0.065 Mq (sl5wl.5r20m3), 0.0065 
Mq (sl5wl.5r20m4). The density structures of the pre-SN 
models are shown in Figure [3] 

Figure [T2] shows the bolometric LCs. The dependence 
on CSM mass in the case a = 1.5 is similar to that of the 
case a = 2. Looking into the bolometric LCs with the same 
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Figure 11. SEDs of the models with different radii at the epoch 
of the maximum bolometric luminosity. The mass-loss rate of the 
models is 10~ 2 Mq yr _1 . 



CSM mass, the LCs with a = 1.5 are fainter until around the 
bolometric peak and then become brighter. This is because 
in the case of a = 1.5, CSM is denser outside and thinner 
inside compared to the case of a — 2 (Figure [3]). Since the 
kinetic energy is more efficiently converted to the radiation 
energy with denser CSM, LCs from shallower CSM become 
brighter at later epochs. Although the luminosity of LCs is 
affected by density slopes, the durations of the round phase 
and the epochs of the sudden drop in the cases of a = 1.5 are 
similar to the cases of a = 2 and are not strongly affected 
by density slopes. 

The multicolor LCs are similar to those shown in Figure 
[7] The NUV absolute Vega magnitudes become as bright as 
-22.8 mag, -21.7 mag, -20.7 mag, and -20.0 mag for the 
models sl5wl.5r20mle3, sl5wl.5r20m2e3, sl5wl.5r20m3e3, 
and sl5wl.5r20m4e3, respectively. 



3.2.5 Dependence on Explosion Energy 

In this section, we look into the effect of the explosion energy 
on LCs. As the higher explosion energy leads to the higher 
kinetic energy of the SN ejecta, the luminosities of the LCs 
during the IPP are expected to be higher with the higher 
explosion energies. The fixed parameters in this section are 
the mass-loss rate (10~ 2 Mq yr" 1 ), the CSM radius (2x 10 15 
cm), the density slope (a = 2), and the progenitor (sl5). The 
explosion energies we adopt are lxlO 51 erg (sl5w2r20m2el), 
3 x 10 51 erg (sl5w2r20m2e3), 5 x 10 51 erg (sl5w2r20m2e5), 
and 7 x 10 51 erg (sl5w2r20m2e7). The density structure of 
the models in this section is the same as that of the model 
Sl5w2r20m2 (Figure [3). 

Figure [13] is the bolometric LCs of the explosions with 
different explosion energies. As expected, the LCs become 
brighter with higher explosion energies. Since the shock 
propagates faster in a higher energy model, the model has a 
shorter rising time and a shorter duration. However, the LCs 
are less sensitive to explosion energies than mass-loss rates 
and radii. This is because the CSM parameters have direct 
effect on diffusion timescales, while the explosion energy de- 
termines the strength of the shock wave where photons are 
emitted. 
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Figure 12. Bolometric LCs of the models with different CSM 
density slopes and mass-loss rates. The bolometric LCs shown 
with solid lines have CSM density slope of p oc r~ 1,5 and those 
shown with dashed lines have CSM density slope of p oc r~ 2 . 'sl5' 
in the model names are omitted in the figure. The bolometric LC 
of the model s!5e3 (no CSM) is also shown for comparison. 
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Figure 13. Bolometric LCs from the model sl5w2r20m2 with 
different explosion energies. 



The multicolor LCs of all the models shown in this sec- 
tion are similar to those of the model sl5w2r20m2e3 in Fig- 
ure [7] with different rising times, durations, and brightness. 
The rising times and durations are the same as those of the 
corresponding bolometric LCs and the NUV absolute Vega 
magnitudes become as bright as —21.0 mag, —22.4 mag, and 
-22.8 mag for the models sl5w2r20m2el, sl5w2r20m2e5, 
and sl5w2r20m2e7, respectively. 



3.2.6 Dependence on 

In the previous sections, the progenitor model inside CSM 
is fixed to sl5. As LCs during the IPP are powered by the 
shock interaction of SN ejecta and CSM, the effect on LCs 
due to the difference in progenitors (RSGs) inside is ex- 
pected to be small. To confirm this, we calculate the LCs 
of the models which have different progenitors but the same 
CSM parameters. The CSM parameters fixed are the mass- 
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Figure 14. Density structures of the pre-SN models with differ- 
ent progenitor models inside. 



Figure 15. Bolometric LCs of the models with different progen- 
itors inside. 



loss rate (10 -2 M Q yr _1 ), the CSM radius (2 x 10 15 cm), 
and the density slope (a = 2). The explosion energy is also 
fixed to 3 x 10 51 erg. We use Type IIP SN progenitors sl3 
(sl3w2r20m2), sl5 (sl5w2r20m2), sl8 (sl8w2r20m2), and 
s20 (s20w2r20m2) to see the effect. Figure [TJ shows the 
density structures of the models. We note that the CSM 
density structures of the models are slightly different from 
each other because the radius where CSM is connected to 
the central progenitor is different depending on the pre-SN 
models. 

Figure [15] shows the LCs with different progenitors. 
Roughly speaking, the LCs are similar to each other because 
all the progenitors are RSGs and properties of progenitors, 
like radii or density structures, are not dramatically different 
from each other. There are slight differences in the maximum 
luminosities and the epochs of the sudden drop in the LCs 
due to the slight difference in those properties of progenitors. 

The multicolor LCs are also similar to those of the 
model Sl5w2r20m2e3 (Figure [7j. The NUV absolute Vega 
magnitudes become as bright as —22.0 mag, —21.8 mag, and 
-21.9 mag for the models sl3w2r20m2e3, sl8w2r20m2e3, 
and s20w2r20m2e3, respectively. 



4 DISCUSSION 

In this section, we compare our LCs obtained in the previous 
section to that of UV-bright Type IIP SN 2 009kf whose LC 
is su ggested to be affected by dense CSM (Botticcll a et al.l 
l201fj) . We show that the LC of SN 2009kf is actually re- 
produced by the LC models with dense CSM and we get 
constraints on the state of the CSM around the progenitor 
of SN 2009kf at the pre-SN stage. 

4.1 Observations of SN 2009kf and Light Curve 
Modeling 

SN 2009kf was discovered by Pan- STAR RS 1 survey and 
obser ved by GALEX using the NUV filter (|Botticella et all 
I2010T I. The observations of SN 2009kf by GALEX revealed 
its distinguishing features: SN 2009kf was continued to be 
bright in NUV for more than 10 days and it was also bright 



in the optical bands during the same period. This feature 
is difficult to be explained by Type IIP SN models without 
dense CSM. This is because, after the shock breakout, the 
UV LCs of Type IIP SNe without dense CSM decline rapidly 
due to the adiabatic cooling of the ejecta and the absorptions 
by iron group elements. Therefore, the optical brightness of 
Type IIP SNe without dense CSM increases as the ejecta 
cools down with the decreasing UV brightness (Figure [2}. 
lUtrobin. Chueai. fc Botticellal (|2010h try to model the LC 
of SN 2009kf without dense CSM and they find that large 
explosion energy (2 x 10 52 erg) is required to obtain the LC 
of SN 2009kf. 

One big uncertainty is in the extinction of SN 2009kf. 
Although t he galactic extinction is negligible (E (B — V) = 
0.009 mag; ISchleeel. Finkbeiner. fc Davis! fl99ct ) . the host 
extinction is estimate d as E(B — V) = 0.32 ± 0.5 mag 
|Botticella et al.l l2010h . The large uncertainty in the host 
extinction makes it difficult to estimate its absolute NUV 
magnitudes. The redsh ift of the host galaxy is 0.182 ± 0.002 
ijBotticella et al.ll2010l ). 

In Figure 1161 we show a comparison of the model 
Sl5w2r20m2e3 with the LCs of SN 2009kf. The SED 
derived at each time step in STELLA is reddened with 
the host extinction assuming the extinction law of 
ICardelli. Clayton, fc Mathlsl (|l989T ) and then shifted to the 
redshift 0.182 by using the standard ACDM cosmology with 
H = 70 km s" 1 Mpc" 1 , D. M = 0.3, and S1 A = 0.7. The host 
extinction we apply here is E(B — V) = 0.22 mag. The left 
panel of Figure [T6l shows the LCs in the first 60 days since the 
explosion in the observer frame. The solid lin es are the model 
LCs a nd the points are the observations by iBotticella et all 
(2010). The epochs of observations are shifted arbitrarily. 
There is little contribution of 56 Ni produced by explosive nu- 
cleosynthesis on LCs at these epochs. The model LCs are in 
good agreement with the observations. Especially, the char- 
acteristic observational feature that the NUV LC and the 
optical LCs are bright at the same epochs is well-reproduced, 
as well as the NUV brightness. In the right panel of Figure 
1161 the LCs at the later epochs are also shown. Since the 
explosive nucleosynthesis (amount of 56 Ni produced), the 
progenitor model inside (mass of hydrogen layer), and the 
degree of mixing of 56 Ni in H-rich layer mainly affect the 
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Figure 16. Comparison of the model LCs from sl5w2r20m2e3 with the observed LCs of SN 2009kf. Each point is observed LCs and 
the solid lines are the model LCs shifted to the observer frame. The epochs of the observation points are shifted arbitrarily. The model 
LCs are obtained by reddening the SEDs obtained by numerical calculations with the host extinction of E(B — V) = 0.22 mag and then 
adopting the redshift 0.182. Magnitudes of NUV are shifted by 3 magnitudes in this figure. 



LCs at the epochs after the IPP (e.g., iKasen fc Wooslev 
20091 ; IDessart. Livne. fc Waldmanl 1201 dl; IDessart fc Hilliei 
20111 ; iBersten. Benvenuto. fc Hamuvl 1201 if ) , modeling this 
part of the LCs is out of our scope. Especially, if we include 
such an effect, photosphere is expected to be located outer 
than our model and photospheric velocity becomes faster. 
The model shown in Figure [16] has the explosion energy of 
3 x 10 51 erg and the evolution of the photospheric velocity 
is shown in Figure [5] with the observed Ha line velocity. Al- 
though our model LCs have slower photospheric velocity, it 
is expected to increase as is mentioned above. Large amount 
of 56 Ni production is expected from the long plateau phase 
after the IPP. In addition, the observational facts that the 
bolometric luminosity at the plateau phase is very high and 
the Ha line velocities of SN 2009kf are very large also indi- 
cate large explosion energy of SN 2009kf. Note that the LCs 
during the IPP are not sensitive to the central progenitor 
model (Section 13. 2. 6[) and progenitor models other than sl5 
can also work. 

4.2 CSM around the Progenitor of SN 2009kf 

Given that the only observational data we are able to com- 
pare with our model LCs are those during the IPP and that 
the difference in progenitors inside (RSGs) does not have 
much effect on LCs during the IPP (Section 13.2.61 Figure 
I15[) . it is difficult to constrain the ZAMS mass of the pro- 
genitor of SN 2009kf. However, as we show in Section 13.21 
the LCs of the IPP is strongly affected by the CSM parame- 
ters and we can get constraints on them. In this section, we 
try to make constraints on the CSM around the progenitor 
of SN 2009kf at the pre-SN stage. 

First, as discussed in Section [3. 2. II the mass-loss rate of 
the progenitor should be larger than ~ 10~ 4 M© yr~* to see 
the effect of the CSM on LCs. The duration of the UV- bright 
phase of SN 2009kf indicates that the CSM radius is larger 
than ~lx 10 15 cm (Section l3~2~3l Figure [9]). Although a 
model with low explosion energy and a small CSM radius can 
result in LCs with a duration similar to that of SN 2009kf, 
the high bolometric luminosity and the large line velocities of 



SN 2009kf are difficult to be reproduced by such low-energy 
explosion. 

The absolute NUV magnitude of the UV-bright phase 
is required to make further constraints on CSM. However, 
due to the uncertainty in the host extinction, it is difficult to 
get the absolute NUV magnitude of SN 2009kf. Therefore, 
what we can confidently conclude is that SN 2009kf has a 
LC naturally explained by the CSM interaction (long-lasting 
UV brightness during the period when it is also bright in 
optical) and the mass-loss rate of its progenitor should be 
higher than 10 -4 Mq yr _1 just before or at the time of 
the explosion to see such effect on the LCs, assuming that 
the CSM velocity is 10 6 cm s" 1 . The CSM radius which is 
larger than ~lx 10 15 cm is inferred from the duration and 
brightness of the NUV LC. If the host extinction is assumed 
to be E(B — V) = 0.22 mag, the multicolor LC model of 
the IPP shown in Figure 1161 which has the CSM with the 
mass-loss rate 10 -2 Mq yr _1 and the radius 2 x 10 15 cm, is 
consistent with the UV-bright phase of SN 2009kf. We note 
that the density slope of the CSM does not have much effect 
on the duration of the IPP as shown in Section [3.2.41 and we 
cannot constrain the density slope with our model. 



4.3 Progenitor of SN 2009kf and Its Extensive 
Mass Loss 

The high bolometric luminosity at the plateau phase and 
very high line velocities of SN 2009kf imply that SN 
2009kf had a high explosion energy. In addition, the long 
plateau phase indicates that the amount of 56 Ni produced 
by the explosive nucleosynthesis is large. Since Type IIP 
SNe from higher ZAMS mass progenitors te nd to be more 
energetic and produce more 56 Ni (e.g., iHamuvl [2003; 
lUtrobin fc Chugai| [2009). the progenitor of SN 2009kf is in- 
dicated to be one of the most massive RSGs and may come 
from the high mass end of RSGs. Therefore, it is indicated 
that the most massive RSGs can have a very high mass-loss 
rate and that such extensive mass loss can occur just before 
the explosions of the massive RSGs. 

Not only LCs but also spectra can be affected by the 
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Figure 17. Left: Possible relation between ZAM S mass and SN typ es. We set the maximum ZAMS mass of Type IIP SNe as M^ p , 
which is observationally suggested to be ~ 17 Mq jSmartt et a,l.l[200gT l . Dense CSM may be left at the time of SN explosions at a ZAMS 
mass range above M™j p . If there is enough CSM, the SNe can be Type Iln. On the other hand, there can be a small ZAMS mass range 
from to M" : p + 5 (with small 8), in which there is not enough CSM for the SN to continue to be Type Iln and the SN can be a hybrid 

of Type IIP and Type Iln. SN 2009kf and SN 1987C are possible candidates for such a hybrid SN. The number fraction of SNe (y-axis) is 
determined by an initial mass function. Right: The same as left but plotted with progenitor stars. The explosion of RSGs without dense 
CSM is expected from ~ 8 Mq to Af™. RSGs explode within dense CSM above M^ p . A yellow supergiant is obse rvationally found 
to be the progenitor o f Type IIL SN 2009kr. The ZAMS mass of the yellow supergiant is estimated as ~ 18 — 24 Mq ( Elias -Rosa etajj 
2010) and 15 +5 4 Mq l lFraser et al.| [2010), which arc consistent with our picture. 



existence of dense CSM. During the IPP phase of SNe from 
RSGs with dense CSM, photosphere is located in CSM and 
very narrow P-Cygni profiles, which are similar to those of 
Type Iln SNe, are expected to be observed (Figure[5]). Then, 
after the IPP, their spectra shift to those of Type IIP SNe. 
In other words, SN 2009kf-like SNe can be a hybrid type 
of Type Iln and Type IIP (Figure 1 1 7f) . If extensive mass 
loss of massive RSGs happen just before their explosions, 
CSM mass and/or radius can be so small that the inter- 
action of dense CSM and SN ejecta ends in early epochs. 
Then, the corresponding SN may be observed as a hybrid 
of Type Iln and Type IIP and SN 2009kf might be classi- 
fied as Type Iln if early spectra were taken. If there is large 
and/or massive enough CSM or a shell exists that is cre- 
ated long before the explosion due to extensive mass loss, 
the SNe may be purely Type Iln and may not show the 
feature of Type IIP SNe (Figure [TT)). Our prediction that 
the early spectra of SN 2009kf-like SNe have narrow lines 
is what clearly differs from theoretic al models suggested by 
lUtrobin. Chugai. fe Botticellal (|2010t ). 

While no spectra of SN 2009kf were taken during the 
IPP, SN 1987C are suggested to have shown such transi- 
tion of the spectra from a Type Iln SN-like blue spectrum 
with narrow hydrogen em ission lines to Type IIP SN spectra 
jSchlegel fc Kirshner|[T99sl ) . The Type Iln SN-like spectrum 
of SN 1987C is taken at 52 days since its discovery when 



SN 1987C could have been approaching to the end of the 
IPP. At 79 days after the discovery, the spectrum showed 
the P-Cygni profile of hydrogen lines and the line velocity 
of Ha was high (6800 km s" 1 ). Although SN 1987C was not 
observed by UV, this transition of the spectra could indicate 
that SN 1987C may be another sample of an explosion of a 
RSG within a dense CSM. The high Ha line velocities indi- 
cates that the progenitor of SN 1987C was a massive RSG. 
Such early observations of SN spectra are important to find 
other candidate SNe which are the hybrid of Type Iln and 
Type IIP. 

Although we show a possible evidence for the exis- 
tence of extensive mass loss of a massive RSG, the mech- 
anism is still uncertain. A candidate for the mechanism 
to ac hieve this kind of mass loss is the pulsations of 
RSG s (|Yoon fc Cantielloll2010l;lHeger et al.lll997l ; lLi fc Gongj 
11994) . As shown by lYoon fc Cantiellol (|2010T L this mecha- 
nism can potentially make a mass-loss rate of a RSG as high 
as 10~ 2 M Q yr _1 . What is more, for the RSGs with ZAMS 
masses around 17 Mq, the pulsations just before the explo- 
sion and the CSM can potentially remain dense until the 
time of the explosion. Coincidentally, th e minimum ZAMS 
mass to cause the pulsations obtained by lYoon fc Cantiellol 
(2010) is roughly the same as the maximum ZAMS mass of 
Type IIP SN proge nitors indicated by observations (17 Mq, 



ISmartt et al.l 120091 ). Such extensive mass loss in massive 



© 2011 RAS, MNRAS 000, HUTS] 



Supernovae from RSGs with Extensive Mass Loss 15 



Type IIP SN progenitors may suppress the upper limit on 
the ZAMS mass of Type IIP SN progenitors because such 
mass loss can take the whole hydrogen layer out of the pro- 
genitors. Mass loss due to nuclear flash may also be a drivin, 
force of the extensive mass loss l|Weaver fc Wooslevl 1 197! 
iDessart. Livne. "fc WaldmarJlioIoh . 

Rareness of UV-bright Type IIP SNe similar to SN 
2009kf can be interpreted as a relatively small ZAMS mass 
range of this event. For example, no Type IIP SNe observed 
with Swift satellite show the long te rm UV-brightenin g 
which is expected by CSM interaction (|Brown et al.l l2009). 
The rareness of UV-bright Type IIP SNe also support that 
they may come from the high mass end of Type IIP SN pro- 
genitors (Figure [17]) . Note, however, that the rareness can 
also be interpreted in different ways. For example, it is pos- 
sible that RSGs generally have extensive mass loss (due to, 
e.g., nuclear flashes) but it usually occurs long before their 
explosions and the mass range of the progenitors which ex- 
perience extensive mass loss just before the explosions is 
small. Future early spectral and UV observations of Type 
IIP SNe are required to get more samples of SN 2009kf-likc 
SNe so that we can make constraint on the mass range of 
SN 2009kf-like SNe and on the driving force of such exten- 
sive mass loss. Observations of nebular phase spectra are 
also i mportant for determining the ZAM S mass of such SNe 
(e.g. , IDessart. Livne. fc WaldmarfeoiOT l . 



5 CONCLUSIONS 

(i) We show that the existence of dense CSM affects LCs 
of explosions of RSGs. This is because (1) shock breakout 
signals are elongated by CSM and (2) SN ejecta is decel- 
erated by CSM. In particular, because of the deceleration, 
kinetic energy of SN ejecta is converted to thermal energy 
which is emitted as radiation and SNe can be brighter than 
usual. The LC becomes bright in UV as well as in opti- 
cal. In addition, the photospheric velocity of early epochs is 
very low because the photosphere is located in CSM at early 
epochs epochs. The most influential parameters of CSM 
are mass-loss rates and radii. The mass-loss rate should be 
higher than ~ 1(T 4 M yr -1 to show the effect of CSM. 
Higher mass- loss rates and/or larger radii lead to longer dif- 
fusion timescales of CSM and thus, longer durations and 
rising times of the LCs powered by the interaction. Density 
slopes and explosion energies also slightly change LCs. The 
difference in Type IIP SN progenitor (RSG) models inside 
CSM are not so sensitive to LCs. (Section [3]) 

(ii) The LCs of Type IIP SN 2009kf, w hich were bright 
in UV as well as in optical in early phases (|Botticella et al.l 
2010), can be explained by the pre-SN models with dense 
CSM. The mass-loss rate of the progenitor of SN 2009kf 
should be higher than 10" 4 M yr" 1 . The CSM radius is 
expected to be larger than ~ 1 x 10 15 cm. The explosion 
energy of SN 2009kf is likely to be very high because of its 
high bolometric luminosity at the plateau phase and its high 
line velocities. The long duration of the plateau phase of SN 
2009kf implies that the large amount of 56 Ni is produced 
by the explosion. 

(iii) The high explosion energy and the large amount of 
56 Ni produced indicate that the progenitor of SN 2009kf is a 
massive RSG. Our results show that massive RSGs are likely 



to experience extensive mass loss exceeding 10~ 4 M yr -1 
just before their explosions. The explosions of such mas- 
sive RSGs with extensive mass loss will be SN 2009kf- 
like SNe. Their spectra show the transition from Type Iln 
SN-like spectr a to Type IIP SN spectra which is observed 
in SN 1987C l|Schlegel fc Kirshnerl Il998h . The existence of 
such SNe indicates that there is some mechanism to en- 
ha nce mass loss of massive RSGs, like pulsations discussed 
bvlYoon fc Cantielld J2010h or nuclear flashes suggested by 
I Weaver fc Wooslevl \l9m . Such a mechanism may reduce 
the maximum ZAMS mass of Type IIP SN progenitors pre- 
dicted by single star evolution modeling (~ 25M ; e.g., 
iHeger et al.ll2003h as low as the observa tionally implicated 
value (~ 17M ; e.g.. ISmartt et alj|2009h (Figure IT?). 

(iv) Future early spectral and UV observations of SNe will 
find other SNe similar to SN 2009kf and provide LCs with 
long time coverage and spectra at the IPP. Large samples 
and detailed observations of SN 2009kf-like SNe can reveal 
the mass-loss mechanism of RSGs as well as the nature of 
SN 2009kf-like SNe. 
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